Two characteristic temperatures for a Bose-Einstein condensate of a finite number of 

particles 
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We consider two characteristic temperatures for a Bose-Einstein condensate, that are related 
to certain properties of the condensate statistics. We calculate them for an ideal gas confined in 
power-law traps and show that they approach the critical temperature in the limit of large number 
of particles. The considered characteristic temperatures can be useful in the studies of Bose-Einstein 
condensates of a finite number of atoms, indicating the point of a phase transition. 
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Experimental achievement of Bose-Einstein conden- 
sation (BEC) in trapped, cold gases of alkali atoms 
0, has stimulated large interest in the physics of this 
phenomenon. Among others, the issue of fluctuations 
in the number of condensed atoms has been a subject 
of intensive theoretical studies, both for noninteracting 
Q0[JSIii0 and interacting gases 0|. For an 

ideal Bose gas, described in the canonical or microcanon- 
ical ensemble, the condensate fluctuations are maximal 
just below the critical temperature. The point of maxi- 
mal fluctuations defines some characteristic temperature, 
which can indicate the occurence of a phase-transition [j| . 
In the systems containing finite number of particles, the 
thermodynamic functions remain analytic for all temper- 
atures, and definition of the critical temperature cannot 
be based on the presence of a singularity in the thermo- 
dynamic functions. In this case, the characteristic tem- 
perature provides an alternative way to localize the point 
of a phase-transition. 

Different characteristic temperature can be defined on 
the basis of the ground-state occupation number. Mean 
ground-state population as a function of the temperature, 
exhibits the presence of the inflexion point close to the 
critical temperature. In this paper we study the proper- 
ties of both characteristic temperatures, demonstrating 
that they tend to the critical temperature in the limit of 
large number of particles. 

The definition of a suitable characteristic temperature 
for finite systems was already considered by several au- 
thors. One of the candidates was a cross-over temper- 
ature. It is the temperature at which the probability 
distribution of the occupation of the condensed state 
changes its character from the one with the maximum 
at zero atoms (above phase transition) to the one with 
the maximum at nonzero (below the phase transition) 
0- fl2l El ■ Another candidate was the maximum in 
the specific heat 0, 0, El . We stress that there are 
other characteristic temperatures possible, all tending to 
the critical one in the limit of N going to infinity and, 
perhaps, one of them may be less difficult to measure in 
the experiment or compute rigorously for the interacting 
system. 



We describe the system of N noninteracting bosons, 
at a temperature T, using the canonical ensemble. The 
canonical partition function Z(N,T), is, by definition 
related to grand canonical one S(z, T), through equa- 
tion: E(z,T) = Y,n=v zNz ( N i T )- Thc S rand canon- 
ical partition function is known in the compact form: 
S(z, T) = ITZot 1 - ze-^-)- 1 where /3 = l/(k B T), and 
e v is the single-particle energy of the level v. Unfortu- 
nately, no such general expression exists for Z(N, T), and 
we calculate Z(N, T) from S(z, T), using Cauchy integral 
formula 



Z(N,T) 



dz S(z,T) 
2ri z N + 1 ' 



(1) 



In the similar manner, it is possible to express the 
ground-state occupation number (No) in the canonical 
ensemble 



(No) = 



1 



dz s(^,r) 



Z(N,T) J 2ni z N + 1 1-z' 



and its fluctuations (S 2 Nq) 



(S 2 N ) = 
(A 2 -7V ) = 



(TV 2 -N ) + {N } -{N } 2 , 
2 f dz S(z, T) ( 



Z{N,T) J 2ni z N + 1 



1 



(2) 



(3) 
(4) 



Standard way of calculating the contour integrals (JTJ , (J2J 
and Q is to use the saddle point (SP) approximation. 
This method exploits the fact that for large number of 
atoms the integrand is sharply peaked around the most 
probable value, and can be approximated by the Gaus- 
sian function. This approximation, however, ceases to 
be valid below the critical temperature, Tc, of the BEC, 
when the saddle point lies close to the ground-state sin- 
gularity. In this case SP method predicts improper re- 
sults for the condensate fluctuations at low temperatures 
11 S El Alternatively, the Maxwell Demon (MD) en- 
semble [3 can be used, in which the moments of the 
condensate statistics can be easily elaborated from the 
grand canonical partition function of excited subsystem. 
This approach, however, is not suitable for studying the 
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condensate statistics in the transition region, and in the 
case of fluctuations, it does not predict the occurence of 
the maximum 00, El- 

In this paper we apply a variant of the SP method , de - 
veloped quite recently by Holthaus and Kalinovski 10] . 
Following the idea of Dingle , the authors of 0] ex ~ 
elude the ground-state term from the Taylor expansion 
around the SP and perform further calculation with the 
special treatment of the ground-state term. In this way, 
they calculate mean number of condensed atoms and its 
fluctuations, in the whole range of temperatures, includ- 
ing the most interesting regime of the phase transition. 
According to we rewrite Eq. 10 in the form 



Z(N,T) = 



dz exp(-FQ,r)) 
2m 1 — z 



(5) 



) is 



where F(z, T) = {N + 1) In z + Y^Li M 1 
a tempered function, with excluded contribution of the 
ground-state term. Function F(z, T) should be expanded 
around the saddle point Zo, which is calculated for the 
whole integrand. It fulfills the following equation 



00 -Be,, 

z e 



1 - Zo 



z e 



-8e„ 



N + l. 



(6) 



Expansion of F{z, T) around z$ up to the second-order 
terms, leads to an integral involving a product of the 
Gaussian function with the singular ground-state term 
1/(1 — z). The result of integration may be expressed 
in terms of the parabolic cylinder function [Tfjj ]. In the 
regime of condensation, it can be approximated by rela- 
tively simple expression 



Z(N,P) = e 



_ -F(z Q .T)-l 



(7) 



Eqs. lfTfj|) - ((TT1) . with the saddle points z , z\ and z 2 cal- 
culated from Eqs. iJBJ, iJSJ, and © respectively, deter- 
mine the mean ground-state occupation number (No) and 
the fluctuations (6 2 Nq) for temperatures T < Tc- Com- 
parison with the exact numerical results calculated for 
moderate-size systems containing 10 2 - 10 4 atoms demon- 
strate, that the modified SP method predicts very accu- 
rate values of (No) and (S 2 Nq) . The accuracy of the SP 
approximation grows with the number of atoms, which 
reflects the fact, that for larger N the integrand becomes 
more sharply peaked in the SP. 

Derivation of the analytical results for the considered 
characteristic temperatures, requires to solve, at least in 
some approximation, the implicit equations determining 
the saddle points. The approximate solutions of Eqs. 
©, © and (0 may be found, by expanding their left 
hand side around z = 1, without altering the ground- 
state term z/(l — z). Applying this procedure to Eq. 
©, we obtain 



zo 



1 - Zo 



+ n e (T) + d e {T)(z - 1) ~ N + 1, 



where 



v=0 



and 



-8e„ 



v=0 



-8e v 



-8s„ 



3e„ \l_ e -8e 



(12) 



(13) 



(14) 



After expansion, the SP Eq. 1)12(1 becomes quadratic in 
z . Its solution reads 



The same procedure may be repeated for the contour 
integrals (0 and (0, describing (N ) and (5 2 N ), re- 
spectively. The SP of the former integral, z\ , fulfills the 
following equation 



00 -8e„ 

1, + ]T I ^ = N 



1- Zl 



z\e 



(8) 



while the SP of the latter integral, z\, has to be deter- 
mined from 



Z-2 



1-Z 2 



St 

i/=i 



z-ie 



-8e„ 



z 2 e 



-8c 



= N-1. 



(9) 



Similar procedure, as that used for derivation of J7J), ap- 
plied to the integrals 0) and (0), yields the following 
results 0] 

(JV ) - 2-5^e*'<*» T >- F <* 1 ' r >- 1 ) (10) 
1 - Zl 

(N 2 -N a ) = 9(y^7) e F(zoJVF(z2 ' T) ~ 2 . (11) 



zo 



no 



2d, 



y/n 2 + 44 



2d, 



(15) 



where rio(T) — N — n,(T) + 2. The same procedure 
repeated for the two remaining SP Eqs. (0 and 10, 
leads to the following solutions 



Zl 



Z2 = 



no + 2d, — y/n Q l + 8d e 
2d, 



no + 2d, — y/riQ + 12d, 
2d, 



(16) 
(17) 



The distance of the saddle points z\ and z 2 to zo, in the 
complex plane, is very small, which may be easily seen 
from Eqs. (Fft, OH, (TJJ. Hence, F(zi,T) and F(z 2 ,T), 
appearing in Eqs. I|10|l and l|ll|) . may be calculated by 
means of the Taylor expansion around z — Zq- This will 
allow to further simplify Eqs. 1|10|) and 1)110. In order to 
preserve sufficient accuracy required for calculation of the 
characteristic temperature, we perform expansion up to 
the second order terms in the small parameters {z\ — zq) 
and (Z2 — zq) for F(z\,T) and F(z2,T), respectively. In 
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FIG. 1: Canonical fluctuations of the condensate occupation 
number for TV = 10 4 atoms. Thick solid line depicts the exact 
numerical result. Thin solid line corresponds to the numerical 
result calculated from Eqs. HUH and 11 U derived in the SP 
approximation. Dashed lines represents the analytical result 




the case of F(z\, T), we obtain 



F( Zl ,T) ~ F(z ,T) 



2 ( d *( T ) 



Zl - Zp 
1-Zq 



1 



ZQ 



Z2 ~ Zq 

zo 



(18) 



where we have applied the following approximation 
f^r(2o) - -d e (T) - Zq/(1 - z ). The small error that 
results from this approximation is, in fact, a higher order 
term, which may be safely neglected. In the same way we 
expand F(z2,T), and substitute the expanded functions 
into (|10|) and i|ll|) . obtaining 



(N ) = 2-^_ x 

x e (^) 2 [T^(l + ^l)+^(T)]-l ; 
2 



(N$-N ) = 9 



Z-2 



1-Z2 



x ev *° 



x 



(19) 



(20) 



We have performed several numerical tests, which proved 
that Eqs. JTDJ) - , with the saddle points given by JTSJl- 
(|T7)l . predict very accurate values of (No) and (S 2 N ). In 
comparison to the initial Eqs. I|ll)[) and (|11|) . they do 
not require numerical solving of equations for the sad- 
dle points. Nevertheless, developed set of equations is 
still too complicated, for the analytical calculation of the 
characteristic temperature. To this end, we expand Eqs. 
((To)l - ljT7)) and i|19 fl -l(2D )l in the power series in small pa- 
rameter iV -1 , assuming that n e ~ N and d e ~ N. After 
straightforward, but tedious calculations we arrive at 



(No) = n Q + -—--, 
2 no 2 



(6 2 N ) = d e 



3 d e 

— - - 2, 
2 no 



(21) 
(22) 



FIG. 2: Temperature of the maximal fluctuations T m (a) and 
the temperature of the inflexion point Tj (b) as a function of 
the number of atoms N, calculated in the canonical ensemble 
for the case of 3D harmonic trap. Thick solid line on both 
plots depicts the critical temperature calculated in thermo- 
dynamic limit. Figure (a): T m given by Eq. (12711 (dotted 
line), is compared with the analytical result accounting also 
for lower order terms (thin solid line), and with exact data 
of numerical calculations (squares). Figure (b): Ti given by 
Eq. 1281 (dotted line) is compared with the analytical result 
accounting also for lower order terms (thin solid line), and 
with exact data of numerical calculations (triangles). 



where in the final result we include the terms of the order 
of O(l), and higher. In the thermodynamic limit, only 
the leading order terms are important: no in Eq. 121|) 
and d e in Eq. 122[) . Those terms can be derived from the 
MD ensemble The inclusion of the remaining terms 
is, however, crucial for the occurrence of the inflexion 
point in (No), and the maximum in (S 2 Nq). 

Figure^Jshows the fluctuations (S 2 No) for an ideal gas 
of N — 10 4 atoms, confined in 3D harmonic trap. The 
exact numerical result calculated by means of recurrence 
relations, is compared with predictions of Eqs. (|11|) and 
()22p. In the case of Eq. I|llf) . the saddle points Zq, z\, and 
Z2, are calculated numerically from Eqs. ©, and 
respectively. It is clearly seen that Eq. I jllj l obtained in 
the SP approximation, fits very well the exact numerical 
curve, apart from the region close to the maximum, where 
a small discrepancy is present. This difference, however, 
does not influence significantly the temperature of the 
maximal fluctuations. On the other hand, in a wide range 
of temperatures Eas. (|ll|) and \\1'2\\ gives almost identical 
values of (5 2 No). Their predictions are different in the 
vicinity of the critical temperature, where no — and 
Ea. (|22|) becomes invalid. It should be stressed, however, 
that this wrong behavior at Tc does not influence the 
position of maximum, which remains exactly the same 
as in ltTT|) . 

Now, we turn to the calculation of the temperature 
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of maximal fluctuations T m and the temperature of in- 
flexion point Ti. We consider an ideal Bose gas confined 
in D-dimensional trap with power-law single-particle en- 
ergy spectra. The condensate fluctuations in such sys- 
tems have been studied Q, 0; in this paper we adopt 
the notation of 8]. The energy spectra is given by 
£{„.} = A^ i=1 ff, where A is the energy gap between 
the ground-state and the first excited state, a is an expo- 
nent depending on the shape of potential, and z/j are the 
integer numbers counting the excitations in i-th carte- 
sian direction. Here, we assume the case of spherically 
symmetric potential, stressing that our analysis can be 
extended for general case of an anisotropic trap. More- 
over, we are interested in the case D > a, in which the 
condensation occurs below the critical temperature Th 
given by k B T c /A = (N/C(d))^ d T(l + l/a)-° |f, 
where d = D/a, and denotes Rlemann's Zeta func- 
tion. In power-law traps, functions n e (T) and d e (T) are 
given by Hill 



"AT) = V ( J 



(23) 
(24) 



where we retain only the leading order terms. In Eq. 
(|24|l . the dimcnsionlcss coefficients 7 and a are given by 
a = ((d — l)/£(d), 7 = d in the case of d > 2, while for 
1 < d < 2 they are given by a = T(l + l/a) D - 2a ((3 - 
d)r(cfr 1 C«T 2/<i , 7 = 2. In the border case: d = 2, 
d e (T) exhibits logarithmic dependence on T, however, 
we will not consider this special case. We substitute Eqs. 
J2HJ) and (J2U into Eqs. and and determine the 



characteristic temperatures from: 
0, obtaining 



dT 



d 2 {N ) 

-dT2-\T=T, 



Ti. 

T 



1 



(3a) 



1/3 



( T /d)-2 

N — 3 — . 



j \ 1/3 

ad \ ^ (■ < /d)-2 
N — 3 — . 



d \d- 1 



and 

(25) 
(26) 



is exactly the same for T m and Ti. In the considered 
regime of parameters "f/d < 2, and both characteristic 
temperatures approach Tq when N — > oo. In the specific 
case, of 3D harmonic potential (d = 3, 7 = 3), Eqs. i|25|) 
and H26fl take form 



T 

Ti 



i-(^L) 1/3 n-^ 



(27) 
(28) 



It is worth stressing that the power-law depence on N 



Figure |21 shows the dependence of T m (upper plot) and 
Ti (lower plot) on the number of atoms N, in the case of 
3D harmonic trap. It compares the predictions of Eqs. 
(|25|l and (|26() with exact results of numerical calculations. 
In addition it also presents the analytical result, derived 
with the inclusion of the lower order terms in n e (T) and 
d e (T). In this case we obtain the same iV -1 / 3 dependence 
of T m and Ti as in Eqs. (EH and (j2Hl, but with a dif- 
ferent prefactors equal to C(2)/2C(3) 2 / 3 + (C(2)/9C(3)) 1 / 3 
and C(2)/2C(3) 2 / 3 + (C(2)/18C(3)) 1 / 3 , respectively. From 
Fig. we see that Eqs. (H7J) and J2BI) predicts slightly 
different values in comparison to the exact data. Nev- 
ertheless, the analytical curve accountinng also for lower 
order terms fits very well the numerical results. 

In conclusion, we have studied two characteristic tem- 
peratures for a Bose-Einstein condensate, defined by the 
point of maximal fluctuations in the ground-state occu- 
pation number and by the inflexion point of the ground- 
state occupation number. For a wide class of power-law 
traps, we have calculated both characteristic tempera- 
tures, showing that they approach the critical tempera- 
ture in the limit of large number of particles. The numeri- 
cal calculations performed for 3D harmonic trap, reveal a 
good agreement with our analytical results. The concept 
of the characteristic temperature may be useful in the 
studies of the BEC with a finite number of atoms, where 
the characteristic temperature can indicate the occurence 
of a phase transition. 
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